arXiv:1505.05074v2 [cond-mat.mes-hall] 10 Aug 2015 


Appearance of effective surface conductivity — an experimental and analytic study 

Jakub Lis0 Mateusz Wojtaszek, Rafal Zuzak, Bartosz Such, and Marek Szymonski 
Center for Nanometer-Scale Science and Advanced Materials (NANOSAM), 

Faculty of Physics, Astronomy and Applied Computer Science, 

Jagiellonian University, ul. St. Lojasiewicza 11, 30-348 Krakow, Poland 
(Dated: August 11, 2015) 

Surface conductance measurements on p-type doped germanium show a small but systematic 
change to the surface conductivity at different length scales. This effect is independent of the 
structure of the surface states. We interpret this phenomenon as a manifestation of conductivity 
changes beneath the surface. This hypothesis is confirmed by an analysis of the classical current 
flow equation. We derive an integral formula for calculating of the effective surface conductivity as a 
function of the distance from a point source. Furthermore we derive asymptotic values of the surface 
conductivity at small and large distances. The actual surface conductivity can only be sampled close 
to the current source. At large distances, the conductivity measured on the surface corresponds to 
the bulk value. 


I. INTRODUCTION 

Surface conductance measurements at micron and sub¬ 
micron scales belong to popular experimental routines 
characterizing semiconductor surfaces, see Ref. [M The 
motivation comes from microelectronics @ as well as 
from molecular electronics projects, and the results de¬ 
pend mainly on the system under investigation. The un¬ 
derlying physics may be determined mostly by surface 
states [j, I 3 1 by some bulk features or a mixture of these 
two limiting cases In the first case, quantum trans¬ 

port theo^ may be applied to understand the results, 
e.g. Ref. y. While in the latter case, classical conduc¬ 
tance theory is effective. Surprisingly, the understanding 
of this well-established theory in the context of surface 
measurements is poor, making even a qualitative inter¬ 
pretation of experimental data rather difficult. 

Conductance measurements on beveled surfaces have 
been used successfully to reconstruct the doping profile 
below the surface with a spreading resistance analysis [1- 
These methods, however, rely on refined numerical 
algorithms and offer no analytic insight. As such, they 
are an important tool in quantitative analysis but poorly 
contribute to our physical comprehension. 

Recently, we have reported surface conductance mea¬ 
surements obtained for germanium samples with an 
atomically clean (001) surface m- The results were 
interpreted within the classical theory without any ref¬ 
erence to the surface states. We observed a slight but 
systematic change of the measured conductivity when 
varying the distance between current sources. Our as¬ 
sumption was, that this was due to the change of the 
conductivity beneath the surface as a result of Fermi level 
pinning, a phenomenon resulting in variations of charge 
carrier densities near the surface. Such modeling of the 
subsurface region was numerically explored for the sur¬ 
face measurements on Si(lll)[|. As we could not find 
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any reference discussing changes of the surface conductiv¬ 
ity at different length scales, we decided to address the 
problem in more detail. Surprisingly, general analytic 
conclusions can be made based on a concise formula for 
the electrostatic potential profile due to the current point 
source. We show that a conductivity change in the bulk 
can not only make surface conductivity appear distance- 
dependent, but it can even vary the asymptotic behavior 
of the electrostatic potential. Furthermore, the deviation 
from the bulk conductivity becomes experimentally inac¬ 
cessible to surface measurements at large length scales. 

Our paper is organized as follows. In sect. |n] we review 
the experimental data motivating our considerations. We 
reanalyze them in Appendix [C] in light of our analytic 
developments. A brief discussion of the classical current 
flow equation is given in sect. IIII Al Next, in sect. ME 
we outline the so-called band bending phenomenon that 
motivates our model of the conductivity varying with the 
distance from the surface. In sect. m we outline our an¬ 
alytic findings while more rigorous treatment and tech¬ 
nical details are discussed in Appendices |A] and |B] . We 
also demonstrate numerical results for two simple cases. 
Finally, in sect. El we comment on the possible finite 
size effects of the physical electrodes. We conclude with 
a summary and we point out several interesting issues. 


II. EXPERIMENTAL DATA 

We report surface conductivity obtained from surface 
conductance measurements for two p-type Ga-doped ger¬ 
manium samples cut from one wafer (MTI corporation, 
cr = 2 — 10 (ricrn)”^). Both the preparation and mea¬ 
surements were done in an ultra-high vacuum system [HI- 
The data are obtained for a well-reconstructed Ge(OOl) 
surface, a hydrogen terminated Ge(OOl) surface, and a 
partially dehydrogenated Ge(001):H surface. Ge samples 
were prepared using the procedure described in Ref. I 111 , 
where some data for the bare surface have also been re¬ 
ported. The hydrogen termination of the surface was 
done following the procedure reported in Ref. [T^. The 
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Figure 1. Conductivity cr as a function of the distance D 
between the current source and drain for three different ger¬ 
manium surfaces. The measurements were done using the four 
point-probe technique in a collinear arrangement outlined in 
details in Ref. P. As some points coalesce, the data are also 
given in Tab. |T] 


partial dehydrogenation was achieved by irradiation of 
Ge(001):H with an electron beam from the Scanning 
Electron Microscope (SEM). As a result, in addition to 
an increased population of unsaturated dangling bonds 
a number of carbon atoms were adsorbed at the surface, 
as verified by Scanning Auger Spectroscopy. 

The investigated surfaces have opposite transport 
properties due to the different electronic structures. 
There are well defined surface bands for the clean 
Ge(OOl) surface [p close to the Fermi level. The 
Ge(001):H surface is isolatiM due to band positions far 
from the Fermi level [l^, H] ■ The surface treated with 
the electron beam is supposed to be a disordered system, 
hence it is weakly conducting at best. Nevertheless, the 
conductivity in all three cases is nearly the same. This is 
an interesting property of the Ge samples that we shall 
address in a future publication. For our considerations 
here, it is important that the surface contribution to the 
current flow can be disregarded. As such, we can use the 
classical description of the current suited for transport 
phenomena in bulk. 

Surface conductance measurements were done using a 
four point probe technique as outlined in Ref. m For 
a given distance D between the current source and drain 
we measured the voltage drop between two additional 
probes. All four electrodes were collinear with the inner 
electrodes probing the voltage drop. The geometry was 
symmetric with respect to the point equidistant to the 
electrodes supplementing the current. For every D the 
voltage drop was measured for several positions of the 
inner electrodes. This allowed us to calculate the con¬ 
ductivity for every D separately using a simple fitting 
procedure. The data were consistent with the model for 
three-dimensional current flow. For n-type doped sam¬ 
ples two dimensional currents were observed, which we 
do not address here neither experimentally nor theoreti¬ 


cally. 

The data shown in Fig. 1 and Tab. |I] were obtained from 
two samples with different tips supplementing the cur¬ 
rent. What is of interest here is the fact that the conduc¬ 
tivity resulting from these measurements changes with 
the distance D. These changes seem to follow a deter¬ 
ministic pattern rather than being of a stochastic origin. 
This led us to believe that these changes are not an arti¬ 
fact but a physical phenomenon. Below, we shall put this 
statement on firm ground. We revisit the experimental 
data in Appendix [Cl 


III. GENERAL FRAMEWORK 


A. Classical model of the current flow 

The classical equation governing the current flow fol¬ 
lows from the Maxwell equations [l5l |. It can be derived 
in three steps. First, as a stationary solution is addressed, 
all time derivatives in the Maxwell equations are set to 
zero. In turn, the electric field E is a curl-free field 

V X E = 0, 

giving rise to the notion of the electrostatic potential 
d)(x), such that electric field is given by its gradient, 
E = Vd). Second, a phenomenological input to the the¬ 
ory is needed to relate the electrostatic field and the cur¬ 
rent density j(x), where a linear relation is usually as¬ 
sumed 


j(x) = cr(x)Vd>(x), (1) 

where the conductivity tr is a scalar. The coefficient a is 
a material property, independent of measurement or the 
sample geometry. Third, from the continuity equation 
(V • j = 0) we obtain the relation 

V (ctV$) = 0. (2) 

This equation is valid in regions where there are no cur¬ 
rent sources or drains, which can be taken into account 
by modifying the right-hand side of the above equa¬ 
tion d, [ 3 . 


D 

Sample A 

Sample B 

[//m] 

Ge(001} 

Ge(001):H 

SEM irr. 

Ge(OOl) 

Ge(001):H 

2 

3.56 ± 0.19 

3.03 ± 0.16 

2.80 ± 0.06 

1.81 ± 0.06 

2.10 ± 0.04 

4 

1.98 ± 0.03 

1.67 ± 0.04 

1.95 ± 0.03 

1.67 ± 0.03 

1.58 ± 0.01 

8 

1.62 ± 0.02 

1.38 ± 0.01 

1.63 ± 0.02 

1.37 ± 0.01 

1.30 ± 0.01 

16 

1.44 ± 0.01 

1.42 ± 0.01 

1.47 ± 0.01 

1.27 ± 0.01 

1.27 ± 0.01 


Table I. Experimental values of conductivity (in (flcm)“^) 
for various surfaces and distances shown in Fig. 1. SEM 
irr. denotes the surface irradiated by the electron beam from 
SEM. 
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B. Conductivity 

The proportionality coefficient cr between the current 
density and the electrostatic field has its microscopic in¬ 
terpretation. For semiconductors it is expressed in terms 
of the material properties, 

a = UeVe + rihVh, (3) 

where ng/^ stands for the electron/hole density and Ve/h 
is their mobility, respectively; see Ref. [i3- 

Band bending is a phenomenon that calls for consid¬ 
ering models with the conductivity a = a{z) changing 
with the distance z > 0 from the surface (at z = 0). 
It is explained in many textbooks, e.g. Ref. [l^, so we 
confine ourselves to a succinct resume in terms of a self- 
consistent calculation scheme, and we illustrate the main 
facts about it in Fig. 2. The phenomenon is due to surface 
states whose energy is determined in relation to bands, 
they may appear below the valence band, above the con¬ 
ductance band, or within the band gap. Assume there 
is no shift of the bands against the Fermi level at the 
surface. As such, the surface states accumulate a large 
electric charge. This gives rise to the creation of an elec¬ 
trostatic field that may be seen in the semiclassical (en¬ 
velope) approximation of a shift of the band positions 
with respect to the Fermi level. But such shifts change 
the electric charge at the surface. In real crystals, we ob¬ 
serve the surface Fermi level at the position arriving as 
a self-consistent solution of the problem sketched above. 
Thus, in equilibrium the semiconductor surface is slightly 
charged and is accompanied by an electrostatic field. Due 
to this field, the relative position of bands with respect 
to the Fermi level varies near the surface, and in this re¬ 
gion the number of current carriers may be different from 
their bulk values. As indicated by eq. m, a variation of 
the conductivity follows naturally. A similar model can 
also be found in Ref. d. More quantitative considerations 
allowing for band bending and carrier density profiles can 
be found in Ref. [H and[2l|. We note, that treating band 
bending and current flow as independent phenomena is 
not entirely correct as both of these phenomena refer to 
the electrostatic potential. However, at sufficiently large 
distances from the current source, the current-associated 
field is so weak that it cannot substantially change the 
band bending. It is also not certain, if the above model 
is straightforwardly applicable to germanium. However, 
strong Fermi level pinning in germanium at the valence 
band is beyond any doubt 0. 


IV. SURFACE CURRENT FROM POINT 
SOURCE 


subsurface layer bulk 



Figure 2. A schematic picture explaining the band bending 
and associated change in the number of electrons and holes. 
The acronym s. s. stands for surface states, rie and Uh for 
the electron and hole density, respectively. The width of the 
subsurface region depends on band bending and bulk doping, 
usually it is in the range given in the picture. The depicted 
situation with upwards band bending corresponds to germa¬ 
nium, where the valence band appears very close to the Fermi 
level. 


Cartesian coordinates x and y and the third one con¬ 
strained to the upper semi-axis z > 0. As already men¬ 
tioned, to take into account current sources and drains 
we need to modify eq. ([2) by adding a source term. A 
common choice is the point-source. As such, we concen¬ 
trate on the following equation, 

V(crVd>) =(5(x-x'). (4) 

To fully formulate the model, we state the boundary con¬ 
dition 


dz 


= 0 , 

z=0 


(5) 


which ensures that there is no current flowing through 
the plane z = 0. Eq. (|4]) is a linear equation and may be 
transformed to a familiar form of a Schrddinger operator 
if we switch from the function $ to a function ^ defined 
by the relation 



The resulting form of eq. o yields 




-A- 




C = ^(x-x'). 


( 6 ) 

(7) 


A. Transformation of the current equation 

Now, we turn to eq. ©• We are interested in solving 
the equation in semi-space, i.e., with unrestricted two 


Further, in the case of interest a depends only on the z 
coordinate. We shall use the representation of the Dirac 
delta in terms of eigenfunctions and eignenvalues of the 
operator L. The eigenproblem factorizes, and in two di¬ 
mensions (x and y) it is trivial. All complexity is in the 
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one-dimensional equation 


I 

dz"^ 


^{k;z) 


k'^ip{k; z). 


( 8 ) 


Now we study this equation, as both its eigenfunctions 
and eigenvalues appear in the solution of eq. 0 and ©■ 
In analogy to the Schrodinger equation, we shall refer to 
the term 


F(z) 




(9) 


as the potential. If the potential vanishes, the corre¬ 
sponding eq. ([HI) is called potential-free. Noteworthy, the 
correspondence between the current equation and the 
time-independent Schrodinger equation is a formal fea¬ 
ture, which will facilitate our reasoning. However, one 
should bear in mind, that the current flow equation we 
deal with is a classical theory of the electrostatic field, as 
shown in sect. mm 

Transformation mixes the electrostatic potential $ 
and the function a. The boundary condition ([5|) then 
reads 


dz'4’{k-,0) = 0 for dzcr{0) = 0, 

dztp{k-0) dza{0) o _/n 

where dzip{k-,0) is dil){k,z)/dz taken at z = 0. 


( 10 ) 

( 11 ) 


B. The potential 

Eq. (|S]) binds the potential and the conductivity ct in a 
highly non-trivial way. To shed some light on its physical 
meaning we introduce a function rj 



Figure 3. Gaussian conductance profile describing a tiny en¬ 
hancement of the conductance at the surface (red dotted line) 
and related potential V (black solid line). The conductance 
profile reads a{z) = (j( 0) [0.833 + 0.267exp(—0.0002z^)]. 








The function reflects how the conductivity, and hence 
the carrier density as implied by eq. 0, changes on the 
logarithmic scale. The potential can be rewritten in the 
following form 

V{z) = idzvf+d^,r,. (13) 

We observe, that the potential is independent of the abso¬ 
lute value of (7, i.e. multiplication of a with any number 
does not alter the potential. As the potential contains 
terms with the first and second derivative of ct, the range 
of the change is important. A given potential value may 
be obtained for a small but abrupt change of ct as well 
as for a large deviation spread on a long distance. If 
the variation of ct is slow enough, the resulting poten¬ 
tial is negligible. We assume it is the case for large 2 
as the bulk conductivity is well defined. The presence 
of the term d^rj appears to be slightly inconvenient from 


Figure 4. Gaussian conductance profile describing a small 
deficiency of the conductance at the surface (red dotted line) 
and related potential V (black solid line). The conductance 
prohle reads ct(z) = (j(0) [l.25 — 0.25exp(—0.0002z^)]. 


a physical standpoint. It requires accurate sampling of 
CT if it is to define a physical model. In this paper, we 
are interested in potentials that could be termed as small 
perturbations around the zero valued potential. 

Fig. 3 and 4 show two potentials corresponding to two 
conductance profiles illustrating the formation of accu¬ 
mulation and depletion layers near the surface. As will 
be shown, physically relevant results are obtained if con¬ 
dition m is assumed. This is, however, at odds with 
the classical [I^ and quantum [T^ modeling of the car¬ 
rier density. As such, these conductivity profiles do not 
appear as results of a theoretical calculation. A profile 
of ct(z) approximately reproducing experimental data is 
given in Appendix [Cl 
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C. Electrostatic potential at the surface 

Not much can be said in general about solutions $ of 
eq. (HI). In the potential-free case, the Green function 
may be calculated exactly, as shown in the Appendix IbI 
However, we can infer a lot about the surface profile of 
the Green function, 

(j){r) = d>(x, y,z = 0; x\ y', z' = 0), 

where radial coordinate r = — x')^ + {y — y'Y 

troduced and the point source is located at (a:', y',z' = Q). 
As shown in Appendix]^ (^(r) satisfies the following in¬ 
tegral formula 



= du , (14) 

where Kq{u) is the modified Bessel function of the second 
kind of zeroth order. This function explodes logarithmi¬ 
cally at the origin and vanishes exponentially for large 
values of u; see Fig. 5. For the assumed potentials, the 
function 0) is bounded and does not quickly vary. 
As such, the asymptotic behaviors of ^(r) may be de¬ 
duced from the above formula. Here we summarize the 
results elaborated in detail in Appendix]^ and informally 
explained in Fig. 5. At this point we do not care much 
about the multiplicative prefactor of 4'{r). In sect. lIVDl 
we will set its normalization to suit the physical context. 
First, we report the short-range limit (r —>■ 0), for which 
we have the following relation: 

—> - [ du Ko{u)^^{oo-,0) = ^. (15) 

Ztt r-s-o r Jq Zr 

The same result would be obtained if there were a uni¬ 
form conductivity cr(0) in the whole sample. This is rea¬ 
sonable, as close to the source the current explores only 
a thin layer in the surface neighborhood where the con¬ 
ductivity can be considered constant. It is the behavior 
of '0(fc; 0) for large k that is important here. 

There are two asymptotic relations for large r de¬ 
pending on whether condition m or (HU is considered. 
For (fTUl) we can assume, that '0(0; 0) ^ 0, and then 

/tIoi 1 

^^0(r) —> - du Ko{u)tlj^{0-,0) = (16) 

Ztt r^oo r Jq Zr 

where (t(oo) corresponds to the bulk value of the conduc¬ 
tivity, cr(oo) = limz_).oo o'(z). It is a very intriguing fact, 
that the asymptotic values of conductivity are given by 
its surface and bulk values, no matter how complicated 
the conductivity profile beneath the surface is. 

The second boundary condition (HU gives rise to a 
modified behavior at large distances 0(r) ^ r~^ if 
dzcr{0) > 0. The case of negative dzcr{0) is not treated 
here for the bound state it has in the spectrum (for more 
details see Appendix We do not know any experi¬ 
mental results demonstrating such a quick voltage drop 


Figure 5. The black solid curve corresponds to the func¬ 
tion A'o(m) with a logarithmic divergence close to u = 0 and 
vanishing exponentially for u oo. Two other curves cor¬ 
respond to some function ip^{u/r-,0) with r = 0.6 (red and 
dashed) and r = 12 (blue dotted). Due to the quick decrease 
of Kq{u), only the region (0; 3.5) contributes significantly to 
the integral, eq. HI- For small values of r, the function 0^ 
gets compressed and it has the limiting value 1 nearly on the 
whole interval (0,3.5). For large values of r, the function 0^ 
is stretched, so that for sufficiently large r it is nearly con¬ 
stant (equal to 0(0; 0)) on the interval. The compression and 
stretching is responsible for the obtained asymptotic behav¬ 
ior. 


on the surface. This is why it appears to be a hint as 
to how the conductivity should be modeled close to sur¬ 
faces. Noteworthy, numerical algorithms often deal with 
thin slices of constant conductivity 0113 parallel to the 
surface. As such, they assume condition ([S]). 


D. Effective surface conductivity 


The asymptotic description obtained for the condi¬ 
tion HU motivates introducing a function / (r) which 
is defined as 


c^(0) 

c^e//(r) 


2 

TT 



du Kq(u)iIi^ (~’^) ■ 


(17) 


To justify this definition, we write the electrostatic po¬ 
tential $ in an infinitesimal range around the source. 




I 1 

27rcr(0) yjrl + ’ 


(18) 


where / is the current supplemented. In this way we get 
a normalization factor for the whole range of r. As a 
consequence of eq. (1141) one can write the following for 
any r and ri 


r jFo(u)0 "(^,O) 

rfi /o°“ du Ko{u)iIj^ 


( 19 ) 
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V. BEYOND THE POINT-SOURCE 
APPROXIMATION 
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Figure 6. Numerically obtained cr^ffir) for the conductivity 
<y{z) described in Fig. [3] The dotted horizontal line corre¬ 
sponds to the asymptotic value o-{z) with 2 : —1 00 . In the 
inset is shown. 
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Figure 7. Numerically obtained cF^ffir) for the conductivity 
ct(z) described in Fig. 2] The dotted horizontal line corre¬ 
sponds to the asymptotic value of the bulk conductivity. In 
the inset i/)^(fc;0) is shown. 


For small ri both eq. m and (fTSl) hold. They may be 
plugged into the above formula resulting in the relation 


A more realistic model of the contact supplying the 
current to the sample involves an area on the surface, 
where both the function $ and the current cr(0)92$ are 
given. This corresponds to the physical situation in 
which the potential and current at the contact are deter¬ 
mined by the supplying electrodes. A detailed analysis of 
the problem goes beyond the scope of the present article. 
While the experimental data were acquired with differ¬ 
ent tips, yielding different contact areas and shapes, the 
results seem not to depend on them. This is why we will 
be satisfied with only a few general remarks. 

The requirement that for a given area $(r, 0) be equal 
to a predefined function, can be satisfied by making use of 
the obtained Green function. For a qualitative discussion 
one can formulate an analogue of the multipole expansion 
known in electrostatics to arrive at quickly vanishing 
(n > 1) terms. 

To take into account the incoming current profile, fur¬ 
ther developments are needed. We observe that the cur¬ 
rent density j corresponds, upon transformation ([6|), to 
the quantity 


j = 

Zy/ (J 

Assuming boundary condition (1101) we can draw an anal¬ 
ogy to the quantum mechanics and interpret the current 
density at the surface as the momentum of the incoming 
particles. This makes it clear that for different j we can 
get different penetration depths, and hence different so¬ 
lutions. The Green function is built as a weighted sum 
of all related eigenfunctions, the larger the eigenvalue 
(energy) the smaller the weight. Thus, the point-source 
formula is dominated by low-energy features. As we have 
observed, it is the long-ranged behavior (r —>■ 00 ) that is 
governed by small eigenvalues (/c —>• 0). Regions that are 
distant enough should not be affected either by source 
finite-size effects or by the value of the current density. 


VI. CONCLUSIONS 


0(r) 


I 1 
2-Kaeff{r) r 


( 20 ) 


Thus, aeff{r) may be viewed as an effective surface con¬ 
ductivity. It reflects the fact that the measured surface 
conductivity varies with the distance from the source due 
to changes in the conductivity below. 

Fig. 6 and 7 show the effective conductivity obtained 
for the potentials shown in Fig. 3 and 4. The results come 
from numerical integration of eq. (1141) with numerically 
calculated functions tfj{k;z). 


Our experimental data show a systematic change of 
the surface conductivity with the distance between the 
probes. To understand the results we have investigated 
the classical current flow equation. We have developed 
a theoretical scheme making it evident that the change 
of the conductivity beneath the surface impacts the re¬ 
sults obtained on the surface, i.e. the well-known formula 
for the current point-source (j) ~ l/or is to be trans¬ 
formed to the form (j) ~ lj(7f,ff{r)r. Asymptotic values 
of (Te// have been found to be the actual (t(0) value close 
to the source and the bulk value far from the source. 
This scheme eases the numerical effort needed to obtain 
the surface conductivity profile. Furthermore, it delivers 
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a framework to discuss and classify different effects ob¬ 
tained in numerical studies such as that in Ref. We 
have confined our considerations to small deviations of 
the conductance parameter. 

There are several interesting open questions mentioned 
in the text, however three issues are of practical interest. 
First, the inverse problem, i.e. the existence of a scheme 
of surface measurements allowing for reconstruction of 
the conductivity profile. Second, capturing the mecha¬ 
nism behind the dimensional reduction [y, [IJ seems both 
interesting and feasible. Third, apparently the condition 
^cr(O) = 0 corresponds to physical observations. It is 
at odds with the way (envelope) wavefunctions are mod¬ 
eled near the surface [T^ . This contradiction needs to be 
clarihed. 
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Appendix A: Calculation Details 

Here we describe in more detail the calculations behind 
the results outlined in sect. lYl To begin with, we write 
explicitly eq. 0 which is to be solved. 

{dl + dl + dl + V{z))S,=5[r-v'), (Al) 

where the potential is calculated from the conductivity 
function 

Viz) = a-V^{z)^y^^. 

The solution of eq. (EB) defines the Green function for 
the differential operator. There are many methods of 
finding Green functions, from which spectral decomposi¬ 
tion is most appropriate to us. We observe that 

the related eigenproblem is solved via 

z) = (p2 + q2 + e^(p^+‘>v)^ik; z), 

(A2) 

where we follow notation introduced in sect. IIVI Both p 
and q are any real numbers, while k may be considered 
positive. 

The Green function may be expressed as 

Gix,y,z;x',y',z') = 
/.oo .oo .oo e^pG-V)+^'}(v-y')Mk;z)i:ik;z') 

L -■ 

(A3) 


we consider real ■!/), so ■i/' = '*/'■ The eigenmodes have to 
obey a uniform normalization condition, e.g. 

dz t/^ik; z)'tp{k'; z) = 5{k — k'). 

Using relation A5 from Ref. [U 
1 /■°“ r°° 



we can integrate out p and q to arrive at the most im¬ 
portant formula in this paper 


POO 

Gix,y,z;x',y',z') = 2 tt / dkKo (kr)'ipik; z)ipik, z'). 
do 

(A4) 

where r = ^/(x — x')"^ + [y — y')^ and Kq stands for the 
modihed Bessel function of the second kind. To proceed, 
the functions '0 and related eigenvaues are needed. This 
can be done for the potential-free case as shown below 
(see Appendix IbI). 

With no harm to generality we set x' = 0 and y' = 
0. As the point source is to be located at the surface, 
we also set z' = 0. We look for the surface profile of 
the Green function, hence we can also put z = 0. The 
function G(x, y, 0; 0,0, 0) corresponds to the sought-after 
surface profile 0(r). Now we can argue for the asymptotic 
formulas given in sect. IIVI The integral to solve is 

POO 

G= dk Ko{kr)'il}^{k-Q). (A5) 

do 

Substituting u = kr, the integral becomes 

G = - [ du Kq (u) 0^ f —; o') . (A6) 

r Jo V r / 

We note that the Bessel function Kq has a logarithmic di¬ 
vergence at the origin and exponentially decays for large 
arguments. Integrating Kq gives 

TT 

du Kq (it) = -, 

see II.4.23 in Ref.l^ As the function 0(A:; 0) is bounded, 
for any desired accuracy e > 0 one can find a value uo, 
for which 




du Kq (u) 0^ (“’ 


< e. 


(A7) 


As a consequence, we can switch to integration over a 
finite interval. 

First we deal with the case of small values of r. To 
this end, we observe that for sufficiently large k the po¬ 
tential may be considered a small perturbation of the 
eigenstate obtained for the potential-free equation. As 
such, ijj{k,z) ~ coskz and 0(fc, 0) —>■ I. Hence, there is 
an argument too such that for any t > too it holds that 








ip(t,0) Ri 1. As such, we can split integration into two 
parts, 

rC = / du Ko {u) (—]0] + / duKo{u)-l. 

(A8) 

The first term vanishes due to shrinking of the integration 
interval with r and finally we are left with the limiting 
value 

^ = |:- (^9) 

Now, we address large values of r. We assume that the 
limit 


lim ipit] 0) = V'(0; 0) 


is finite and does not vanish. Hence, there is again to 
such that for any t < to the difference tp{t;0) — ipiO^O) 
can be neglected. Then 


rC = 


/Q 


^min{wo ,io^} 


du Kq (u) (0; 0) + 
du Kq (m) ■ 


(AlO) 


Now it is the second term which vanishes due to reduc¬ 
tion of the integration interval. As such we arrive at the 
limiting value 


C = 


TT'ip'^ ( 0 ; 0 ) 


(All) 


It is easy to note, that eq. ([5]) has an zero eigenvalue 
corresponding to the function y/a. This function is re¬ 
lated to the freedom left by eq. which admits shifting 
its solutions by any constant. The second independent 
solution ipoiz) assuming k = 0 reads 


pJo{z) = \/a{z) f 

Jo cr(2/) 


where to stands for the maximal argument, for which 
the linear approximation may be considered exact. As a 
consequence, we obtain the relation 



with a different asymptotic current behavior. The nu¬ 
merical value of the integral above yields 



as can be interfered from formula 11.4.22 in Ref. [2^. The 
case with a negative value of m is more complicated. The 
spectrum then has an isolated eigenfunction ~ exp mz 
with a negative eigenvalue and hence the Green func¬ 
tion becomes more complicated than the functions con¬ 
sidered here. 

It is evident from the numerical examples described in 
sect. IIVI that CTe// changes on a slightly larger scale than 
the potential. To remark on finding solutions of eq. dS]) 
for small eigenvalues k, we denote the range of the po¬ 
tential V with d. It corresponds to the distance from 
which the potential may be neglected. Then, in the basis 
of {cosfcz}fc>o the matrix element Vk,p reads 

1 

^k,p =-z / dz V{z)[cos{k + p)z + cos{k — p)z\. 

2 Jo 

(A15) 

For small values of k and p the integral tends to the value 




k ,p— ^0 J Q 


dz V{z), 


(A16) 


if in the range (0,d) functions cos(-) can be considered 
constant. As such, the vectors k obeying the inequality 


kd < O.I 


It vanishes for z = 0 but has a non-vanishing first deriva¬ 
tive. Following the theorems on continuity of solutions 
of a differential equation with respect to a parameter 
we make use of the solution to arrive at 

V-(0;0) = WA, (A12) 

y cr(oo) 

where the denominator appears due to the proper nor¬ 
malization. 

To complete the discussion, we address the case, where 
tp{t-, 0) ^ t for t —>• 0 and m > 0 (see eq. (1611) 1. In such a 
case we arrive at the following relation 


Emintuo.tor} 

+ j_ du Ko {u)tp'^ (y:’®) ’ 


rtor 

/ du u^Ko (u) -I- 

Jo 


(A13) 


may be considered small. This hints at the length scale 
needed to solve the problem numerically - to explore the 
lowest eigenvalue sector one needs to work on segments 
~ 100c? in length. We used the standard Mathematica 10 
routines to diagonalize the operator L (with higher pre¬ 
cision). The segment lengths that we explore are about 
20 000 nm and 40 000 nm. The lowest eigenvalues turned 
numerically unstable. This is why we excluded several 
lowest eigenvalues. '0^(0, 0) was found by extrapolation 
a function obtained as quadratic interpolation of '0^(fc, 0) 
for several lowest remaining eigenvalues. As shown in 
Fig. 6 and 7, the results are close to the exact values 
cr(0)/(T(oo). A more credible and efficient approach to 
finding ip{k; 0) is to solve eq. dS]) for different values of 
k with the boundary condition on a segment larger 
than the range of the potential V. The normalization re¬ 
quires then, that there is a given amplitude of oscillations 
for all solutions tp{k, z) in the potential-free region. 
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Appendix B: Green functions for the potential-free 
equation 


^/’^(fc;0) vanishes for /c —>■ 0. As such, the following 
asymptotic expansion for large r 


For the potential-free case the eigenfunctions for the 
eigenvalues values yield cos(fcz) if condition (|TTH) is to 
be satisfied and cos {kz + if{k)) if condition (HID holds. 
The phase shift (f is given by the formula 

ip = —siqnim) arccos , , (Bl) 

where m denotes 92 (t( 0 )/ 2 (t( 0 ) and sign{'m) is the sign 
of m. For eigenfunctions cos{kz) we can complete the 
integration in eq. (IA4D to arrive at the following formula 


G[x,y,z]x' ,y' ,z') = 

1 

27r / dkKQ {kr) - [cos k{z — z^) + cos k{z z')] . 
Jo 2 


'0 

Using relation 11.4.14 in Ref. [2^ we obtain 

TT^ f 1 

G{x,y,z;x',y',z') = 


(B2) 


2 I + (z- z'Y 




1 


(B3) 


Cr2 + (z + zOV ’ 


which reduces to the familiar formula r~^ at the surface. 
The term with z + z' appears due to the boundary condi¬ 
tions to ensure that no current escapes through the sur¬ 
face. Similar Green functions are known in electrostatics 
to give rise to image charges. 

We note that eigenfunctions cos kz converge to a con¬ 
stant in the limit fc —>■ 0 , which solves potential-free 
eq. ([S|) with k = 0. This is a demonstration of theorems 
on the continuity of solutions of a differential equations 
with respect to parameters. 


Appendix C: Experimental data revisited 

The insight gained from the above analysis of the cur¬ 
rent flow equation prompts a few comments on the ex¬ 
perimental data presented in sect. m The bulk value 
of the conductivity appears different for the two sam¬ 
ples investigated and reads about 1.4— 1.5 ll'^cm”^ and 
1.2 — 1.3 f2“^cm“^ for sample A and B, respectively. In 
both cases it is below the nominal value. The exper¬ 
imental conductivity decreases with the distance from 
the source, see Fig. 1. The trend is reproduced in Fig. 6 , 
which suggests an enhancement of the conductivity at the 
surface. For the p-type doped samples the surface Fermi 
level is located close to the valence band and hence the 
formation of an accumulation layer consistently explains 
the data. 

In Fig. 6 and 7 we observe, that the first derivative of 


<^eff{r) =cr(oo) (^1 - 0 ) 2 r^ " ') 

is natural. We stop at the first correction term and 
rewrite the above formula using an effective parameter 


cr^ffir) = cr(oo) + ^. (C2) 

One should bear in mind, that the terms neglected in 
eq. (ICID can be important if the term with cti becomes 
relevant. Hence, some caution is needed when attribut¬ 
ing a physical meaning to ai. 

To apply the ex pre ssion (IC2D to the experimental setup 
described in Ref.llJwe put the point current source at the 
point ( 0 , 0 , 0 ) and the drain at (D, 0 , 0 ) and calculate the 
voltage drop between probes located at xi = 0 , 0 ) 

and X 2 = (l?i^,0,0). Here the parameter 0 < s < 1 
corresponds to the parameter x in Ref. [ll|. The result¬ 
ing formula for the measured resistance R{D; s), i.e. the 
voltage drop between Xi and X 2 divided by the current 
supplemented to the sample, is poorly informative. How¬ 
ever, it is interesting to expand it in powers oi 1/D 


R{D;s) = 


1 4cris(s^-I-3) 1 


7 rCT(oo)(l — s^) D 7 rCT^(oo)(l — 


(C3) 

We confine our considerations to the two lowest orders 
of expansion written above. For s close to zero they are 
finite, thus the first term (of the order D~^) dominates. 
At this level, the scaling behavior characteristic for three 
dimensions ll| may be seen: the quantity D ■ R{D, s) 
does not depend on D. In the case of experimental data 
such a universality is clearly observed for D equal to 8 
and 16 /im, see Fig. 8 . For smaller D the second term 
D~^) spoils the universality. The reason is, that for 
s close to unity, both terms on the right-hand side of 
eq. (IC3D are singular. The second one explodes quicker 
and dominates for big enough s, no matter how large D 
is. As the second term in eq. (IC3D is due to the subsur¬ 
face variation of the conductivity, the deviations from the 
universal behavior bear information about the subsurface 
region. 

The formula assuming (|C2D for R{D; s) was used to fit 
the experimental data for both samples, see Tab. n As 
a result, estimated values of (t(oo) vary much less than 
for a naive fit behind numbers presented in Tab. n The 
values of cti are usually of the same order, and no clear 
trend is seen, which is probably due to the higher order 
terms. 

Dimensional analysis allows to identify a length parame¬ 
ter 


S 



R. 
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Figure 8. Experimental resistance multiplied by D measured 
for various s for sample A with passivated surface. For 
s < 0.2 all the data coincide. The data for D = 8 fj,m 
and D = 16/im coincide for nearly all values of s. Solid 
lines show numerically obtained data for the proHle (j{z) = 
(4.5 cosh“^ lio + ^) '’■(oo). 


function to model conductivity changes 


cr(cx)) 


a 

coshf 


+ 1 : 


(C4) 


where d, a and tT(oo) are parameters adjusted to the ex¬ 
perimental data. This function satisfies the boundary 
condition and it decays exponentially except for a 
region close to the surface. As such, it a reasonable model 
describing screening of the surface charge [l^. The pa¬ 
rameters d and d have physical meaning while (t(oo) de¬ 
pends on the way one solves the problem (applied nor¬ 
malization condition). As shown in Fig. |8] simulated data 
for d = 120 nm and a = 4.5 satisfactorily correlate with 
the experimental points. Similar results can be obtained 
for d ranging between 80 and 200 nm and for d between 
3 and 7. Interpretation of these parameters within the 
band bending theory, e.g. the Schottky approximation, 
seems to be misleading as it overestimates the doping and 
underestimates the band-bending. 


a characteristic distance where surface corrections are im¬ 
portant. In our case, it appears about 150-500 nm. This 
is consistent with the profile a{z) we apply to reproduce 
our experimental data, see Fig. [5] We use the following 


D 

Sample A 

Sample B 

[fim] 

Ge(OOl) 

Ge(001):H 

Si8iVl irr. 

Ge(OOl) 

Ge(001):H 

2 

1.87 / 8.64 

1.87 / 8.64 

2.30 / 3.97 

1.50 / 1.03 

1.93 / 2.65 

4 

1.73 / 7.53 

1.73 / 7.53 

2.05 / 4.97 

1.66 / 2.90 

1.55 / 4.34 

8 

1.55 / 5.84 

1.55 / 5.84 

1.75 / 6.59 

1.75 / 3.83 

1.42 / 7.27 

16 

1.47 / 17.1 

1.47 / 17.1 

1.47 / 12.4 

1.19 / 40.8 

1.20 / 22.4 


Table II. cr(oo) / ai, i.e. bulk conductivity cr(oo) (in (Ilcm)”^) 
and (Ti (in 10“^"^f2“^cm) from fitting for various surfaces. The 
uncertenity for a{oo) is not larger than 5% while for ai it 
varies from several percent for larger D up to 20% for D = 
2 /im. SEM irr. denotes the surface irradiated by the electron 
beam from SEM. 
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